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ABSTRACT 



Variations in the pulsation arrival time of five independent pulsation frequen- 
cies of the DB white dwarf EC 20058-5234 individually imitate the effects of reflex 
motion induced by a planet or companion but are inconsistent when considered 
in unison. The pulsation frequencies vary periodically in a 12.9 year cycle and 
undergo secular changes that are inconsistent with simple neutrino plus photon 
cooling models. The magnitude of the periodic and secular variations increase 
with the period of the pulsations, possibly hinting that the corresponding physi- 
cal mechanism is located near the surface of the star. The phase of the periodic 
variations appears coupled to the sign of the secular variations. The standards 
for pulsation timing based detection of planetary companions around pulsating 
white dwarfs, and possibly other variables such as subdwarf B stars, should be 
reevaluated. The physical mechanism responsible for this surprising result may 
involve a redistribution of angular momentum or a magnetic cycle. Additionally, 
variations in a supposed combination frequency are shown to match the sum of 
the variations of the parent frequencies to remarkable precision, an expected but 
unprecedented confirmation of theoretical predictions. 

*Based on observations obtained at the Southern Astrophysical Research (SOAR) tele- 
scope, which is a joint project of the Ministerio da Ciencia, Tecnologia, e Inovagao (MCTI) da 
Republica Federativa do Brasil, the U.S. National Optical Astronomy Observatory (NOAO), 
the University of North Carolina at Chapel Hill (UNC), and Michigan State University 
(MSU). 
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Introduction and Formalism 



The orbital motion of a planet hosting star about the system's mass center, often 
referred to as a "wobble", will cause variations in the time it takes for light from the star 
to reach external observers. If the star provides some predictable repetitive behavior, 
like regular pulsations, the wobble will cause systematic differences between the time the 
behavior is predicted to be observed and when it is actually observed. This phenomenon 



has resn 



ted in the discovery of t he first known 



planets ( IWolszczan fc Frail 



1992 



Thorsett et al. 



exop lanet and the rest of the pulsar 



19931), confirmation and detection of a 



growing number of planets with Kepler transit timing variations (see 



Holman et al. 



and s ubsequent Kepler pubhcations), a planet around a sub-dwarf B star 



2007 ) as well as strong d etection limits on several pulsating white dwarfs (IMuUallv et al. 



2008 



Hermes et al. 



2010 



Silvotti et al 



20101 ) and one tentative white dwarf planet candidate (jMullally et al. 



20091 ). However, changes to the actual frequency of the repetitive behavior will also 
manifest as differences between the expected and actual time the behavior is observed. To 
illustrate this, imagine observations of two distant ticking clocks. Clock one keeps poor 
time. Everyday the period of the ticks changes in such a way that clock one appears a 
second slow at noon and a second fast at midnight. Clock two keeps perfect time but is 
moving in a light second orbit in the plane of the observer with a period of a day. While 
these clocks are fundamentally very different, the two clocks are indistinguishable in terms 
of timing alone. It has been asserted via lex parsimoniae (Cecum's Razor), specifically 
for white dwarf pulsators, that any sinusoidal or "planet like" variations in the timing of 
the pulsations are most likely due to orbital motion as there is no known mechanism that 
would conspire to make the actual pulsation frequency vary in such a manner. This paper 
presents conclusive evidence to the contrary: frequency variations in a white dwarf pulsator 
that can mimic the effect of a planetary companion. 
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The dominant cooling mechanism for the hot pulsating DB white dwarfs is thought to 
be the production of plasmon neutrinos. The frequency of the pulsations are determined 
in part by the temperature, so cooling should cause a slow secular change in frequency. 
This effect has bee n well studied with the far cooler, photon cooling dominated DA white 



dwarf G117-B15A f lKepler et al. 



20051 ). If the rate of frequency change of a hot pulsating 



DB white dwarf was measured to sufficient precision, the neutrino productio n rate could 



20041 ). This paper 



be extracted and used to constrain the Standard Model (iWinget et al.l 
presents evidence that there are other processes that cause large secular changes to multiple 
pulsation frequencies in a white dwarf. 

The diagnostic tool used to measure variations in pulsation arrival time and pulsation 
frequency is called the O — C, or "observed minus calculated" diagram. Depending on 
the implementation, "observed minus calculated" can be a misnomer. In many cases it is 
a plot of successive measurements of the absolute phase of a periodic variation assuming 
constant frequency. The ordinate and abscissa of the O — C diagram are both expressed 
in units of time, but are sometimes normalized by the typical period of the behavior into 
units of cycles (often labeled "epochs" on the abscissa), radians, or degrees. The O — C 
diagram can reveal frequency changes that are difficult to detect by direct measurement. In 
addition to the detection of planets around pulsars and other pulsating stars, it was the 
primary diagnostic tool that helped earn Hulse and Taylor earn the Nobel Pr ize in Physics 



for their indirect detection of gravitational wave emission ( iTaylor et al. 



diagram remains a cr itical tool acro s s time domain astron o my, s ee e.g. 



1979). The O - C 



Oksala et al. 



Barlow et al 



feonl); 



MuUallv et al. 



teoOSf ): 



Hermes et al 



fl2012h 



(120121); 



For pulsating white dwarf stars and other short period variables it is convenient to 
describe O — C as a time varying absolute phase. Given some fixed observed frequency, fobs) 
the instantaneous amplitude of a pulsation can be rewritten in terms of the time varying 
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absolute phase, r(t) 



H{t) = Asm{2nUs{t-rm 



rit) can be mathematically described as (for a complete derivation see 



Dalessio et al. 



2OI2I) 



/,\ , ^fobs, 1 

r(t) = to — t - 



ohs 



ohs 



'7(0*' + ^ 



(2) 



where to is an arbitrary offset depending on choice of t = 0, jobs is the observed 
frequency, bjohs is the difference between the observed and actual, time averaged, Doppler 
shifted frequency, ^{f) is the perturbation to the frequency, and g(t) is the perturbation to 
the distance between the object and observer. Note that 7(0) = g(0) = 0. 

Except in the case of continuous observation where the number of elapsed cycles is 
explicitly known, all points in the O — C diagram suffer from the ambiguity that they can 
be shifted up or down by any integer times the pulsation period. A common assumption is 
that O — C has been sampled often enough that 7(t) and q{t) change slowly and smoothly 
during the observational gaps. Under this assumption, the location of the points on an 
O — C diagram can be constrained if there is a single unambiguous configuration that 
creates a smooth, continuous trend. 



2. Observations 



EC 20058-5234 flKoen et al. 



I995I ) has been identified as a prime candidate for 
measurement of the neutrino production rate in a white dwarf (see Section [T]). This 
motivated many of the over 250, 000 PMT (photo-multiplier tube) and CCD (charge- 
coupled device) photometric measurements spanning 1994-2011. EC 2005 8-5234 was a 



primary target of the Whole Earth Telescope campaign XC0V15 in 1997 (jSuUivan et al 
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20081 ). followed by seasonal high speed photometry from the Im McLellan tele scope at 



Mt. John Observatory and the 36 inch SMARTS telescope at CTIO (see e.g. 



Sullivan 



20051). It was a: 



( ISuUivan et al. 



so a tertiary target during XCOV27 and was observed at Magellan in 2003 



20071 ) . A summary of the observations is shown in Table [H 



Table 1:: Journal of observations and binning scheme. 
Note that only the lightcurves from XC0V15 and 
XCOV27 consist of data from multiple sites. 



Observatory 


Date Range 


Instrument 


Observations 


SAAO 


1994/05/14 - 


1994/06/02 


PMT 


5535 


SAAO 


1994/07/06 - 


1994/07/12 


PMT 


6881 


SAAO 


1994/10/03 - 


1994/10/03 


PMT 


1304 


SAAO 


1994/06/12 - 


1994/06/12 


PMT 


1378 


XC0V15 


1997/07/02 - 


1997/07/11 


PMT 


46289 


Mt. John 


1997/10/05 - 


1997/10/06 


PMT 


3162 


Mt. John 


1998/07/24 - 


1998/07/25 


PMT 


7057 


Mt. John 


1998/08/14 - 


1998/08/14 


CCD 


1153 


Mt. John 


1999/09/09 - 


1999/09/13 


PMT 


1517 


Mt. John 


2000/07/06 - 


2000/07/08 


PMT 


8025 


Mt. John 


2000/09/05 - 


2000/09/05 


PMT 


770 


Mt. John 


2001/03/30 - 


2001/04/01 


PMT 


2101 


Mt. John 


2001/09/21 - 


2001/09/24 


PMT 


7710 


Mt. John 


2002/04/12 - 


2002/04/15 


PMT 


8872 


Mt. John 


2002/08/01 - 


2002/08/06 


PMT 


12672 


Mt. John 


2002/09/05 - 


2002/09/09 


PMT 


3814 



Continued on Next Page. 
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Table 1 - Continued 



Observatory 


Date Range 


Instrument 


Observations 


Magellan 


2003/07/11 - 


2003/07/13 


CCD 


2317 


Magellan 


2003/07/25 - 


2003/07/31 


CCD 


8865 


Mt. John 


2003/08/27 - 


2003/09/02 


PMT 


12216 


Mt. John 


2003/09/22 - 


2003/09/23 


PMT 


4694 


Mt. John 


2003/10/31 - 


2003/10/31 


PMT 


1229 


Mt. John 


2004/04/23 - 


2004/04/25 


PMT 


2474 


Mt. John 


2004/05/16 - 


2004/05/16 


PMT 


1599 


Mt. John 


2004/06/09 - 


2004/06/16 


PMT 


14554 


Mt. John 


2004/07/09 - 


2004/07/14 


PMT 


17240 


Mt. John 


2004/08/07 - 


2004/08/10 


PMT 


5072 


CTIO SMARTS 


2004/08/24- 


2004/08/28 


CCD 


2320 


CTIO SMARTS 


2005/09/21 - 


2005/09/25 


CCD 


979 


CTIO SMARTS 


2006/08/31 - 


2006/09/13 


CCD 


5019 


CTIO SMARTS 


2007/08/17- 


2007/09/06 


CCD 


14117 


Mt. John 


2007/07/16 - 


2007/07/17 


PMT 


6281 


CTIO SMARTS 


2008/08/29 - 


2008/09/01 


CCD 


1736 


Mt. John 


2008/08/30 - 


2008/08/31 


CCD 


3092 


XCOV271 


2009/05/18 - 


2009/05/26 


CCD 


4609 


Mt. John 


2010/07/05 - 


2010/07/11 


CCD 


20201 


CTIO SMARTS 


2011/09/08 - 


2011/09/15 


CCD 


3584 



tXCOV27 included observations from SOAR, SAAO, and the CTIO SMARTS .9m. 
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3. Results and Analysis 



3.1. Data Reduction and Preparation 



The PMT ph otometry obtained at SAAO in 1994 was reduced as described in 



Koen et al 



fll995[ ). while the 1997 XC0V15 photometr y and the Mt. John PMT 



photometry was reduced using the methods outlined in 



Sullivan et al. 



fl2008[ ). All the 



Mt. John PMT phot ometry collect ed between July 1998 and July 2007 employed a three 



channel photometer ( ISuUivan 



20001 ) attached to the one meter McLellan telesc ope. All 



CCD image cali bration and aperture photom etry was performed with Maestro (iDalessio 



20101). WQED flThompson fc MuUallv 



20091 ) was us ed to remove spuriou s points and to 



Eastman et al. 



2010h for the CCD 



perform the barycentric correction (for a review see 
data. The barycentric correction applied to the PMT data were performed using the 



same implementation as used in WQED ( IStumpfl 



19801). The photometric measurements 



were divided into 36 individual lightcurves, one for each observing campaign in Table [T] 
The following analysis was also performed by dividing the points into lightcurves for each 
individual night and observing season. The differences were not significant and do not merit 
further discussion. 



3.2. Testing for Stability 



The frequency and amplitude of all pulsation frequencies with amplitudes above 0.1% 
were extracted from two of the largest data sets, XC0V15 in 1997 and the 2007 CTIO 
SMARTS observations. The results are summarized in Table [21 The same spectrum of 



Sullivan et al. 



20081) 



frequencies appears in both years (see the Fourier spectrum from 
with the exception that pulsation frequencies F, J, and L have amplitudes below the 
noise threshold in 2007 and that pulsation frequency E appears to have changed at a 
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statistically significant level. The systematically higher pulsation amplitudes in 2007 are an 
observational effec t due to two com panions within several arc seconds of EC 20058-5234 
(see the Figure in iKoen et al 



19951 ). The PMT observations included all three objects as 
it was impractical in the typically 2" seeing conditions at Mt. John to separate flux from 
the target only. However the CCD data, especially when combined with the better seeing 
conditions at CTIO and the Magellan site, permitted extraction of just the flux from EC 
20058-5234 using synthetic aperture techniques. This effect makes it non-trivial to search 



for amplitude modulation, which, particularly for RR Lyrae stars, is 



Buchler fc KoUath 



mown to be coupled 



20111 ). However, if 



to phase modulation (the Blazhko effect, see e.g. 
the pulsation amplitudes are intrinsically constant, the ratio of amplitudes from one set 
to another should be constant. The ratio of the pulsation amplitudes from 2007 to the 
pulsation amplitudes in 1997 is shown in the last column of Table HI The ratio is consistent 
with there being no amplitude variation with a reduced value of 2.3 with a statistical 
likelihood around 10%. This is the worst reduced value for the model of constant 
amplitude between any two of the data sets in Table [H This value of reduced is slightly 
higher than expected for a constant amplitude. The removal of atmospheric extinction 
and other observational effects could introduce small changes in observed amplitude so we 
do not find this slight disagreement particularly alarming. Even if the amplitude of these 
pulsations frequencies are intrinsically changing, the changes are relatively small and only 
barely distinguishable between our two best data sets. 



3.3. Extracting Precise Frequencies 

Precise measurements of the pulsation frequencies were made using the "bootstrapping" 
procedure. The process is described as follows. The frequencies above the noise threshold in 
the July 2004 dataset were extracted. A sum of sinusoids (one sinusoid at each frequency) 
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Table 2: Pulsation frequencies with amplitudes above 0.1% as directly measured from 
XC0V15 and the 2007 CTIO SMARTS data sets. The standard error in the measure- 
ment of the amplitude of all frequencies was 0.09 mma (1 mma = 0.1%) in the 2007 CTIO 
SMARTS dataset and 0.2 mma in the XC0V15 dataset. Several frequencies from XC0V15 
were not detected in the 2007 data set but would have been had the amplitude remained 
constant. The systematic increase in amplitude from 1997 to 2007 is an observational arti- 
fact. The transition from PMTs to CCD imaging which allows removal of some of the flux 
of the nearby companions. The ratio of pulsation amplitudes from 2007 to 1997 is consistent 
with having a constant amplitude with a reduced value of 2.3 and a likelihood of around 
10%. 



Label 


Period (s) 


Frequency (uhz) 


Amplitude (mma) 




(approx) 


1997 


2007 


1997 


2007 


Ratio 


A 


134 


7452.18(.05) 


7452.27(.02) 


1.67 


2.2 


1.3(.l) 


B 


195 


5128.53(.03) 


5128.71(.02) 


2.51 


3.3 


1.31(.09) 


C 


204 


4902.14(.05) 


4902.14(.02) 


1.47 


2.6 


1.8(.2) 


D 


205 


4887.84(.03) 


4887.91(.02) 


2.62 


3.2 


1.22(.09) 


E 


257 


3893.142(.009) 


3893.239(.005) 


8.37 


11.3 


1.35(.03) 


F 


275 


3640.37(.07) 




1.13 






G 


281 


3559.032(.009) 


3559.033(.005) 


8.45 


10.5 


1.24(.03) 


H 


287 


3489.07(.05) 


3489.05(.02) 


1.58 


2.2 


1.4(.l) 


I 


333 


2998.67(.03) 


2998.71(.01) 


3.01 


4.3 


1.43(.08) 


i J 


350 


2852.36(.09) 




1.40 






K 


525 


1903.61(.04) 


1903.49(.02) 


1.83 


2.44 


1.3(.l) 


L 


540 


1852.53(.04) 




1.86 
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were fit to the July 2004 dataset using nonlinear least squares and allowing each frequency 
to converge to the optimum value. After the fit converged, additional data were added to 
the set, ensuring that each of the frequencies had been measured to sufficient precision so 
that if the frequency had remained constant, there would be no ambiguity in the number 
of cycles between any two data points in the combined set. This ensures that the linear 
term of equation [2] will not dominate the O — C diagram and that the frequencies are 
clearly separated from the aliases of the combined set. After the new data were added, the 
frequencies were re-fitted. This process was repeated until most of the 2003 and 2004 data 
were combined into a single set. The amplitudes of pulsation frequencies F and L were well 
below the noise threshold and were not considered. The combined set of all data between 
September 2003 and August 2004 gave precise enough frequency measurements to merge 
the entire data set. However, merging the September 2003 - August 2004 data set with 
data from August 2003 or 2005 decreased the precision of our frequency measurements. 
This is an indication that the phase/frequency of the pulsation frequencies are changing 
by a substantial amount over that time span. Further combination was abandoned beyond 
the September 2003 - August 2004 pulsation frequency measurements. The results are 
summarized in Table [31 



3.4. Building the O — C Diagrams 



A sum of sinusoids at the 10 frequencies in Table |3] were linearly fit to each data set. 
The calculated absolute phase was plotted as O — C for each pulsation frequency. The 
O — C diagram of all ten pulsation frequencies showed structure that was not consistent 
with scatter. It was also apparent that for some frequencies, O — C was changing by a 
large fraction of the period, and that none of the O — C diagra ms were consistent w ith the 



parabolic trend of a simple neutrino plus photon cooling model ( iBischoff-Kim 



20081 ) 



It was assumed (see Section [T]) that for each pulsation frequency, 7(t) and q{t) change 
slowly and smoothly during the observational gaps. An unambiguous smooth trend in the 
O — C diagrams of pulsation frequencies A, D, E, G, H, and I was clearly discernible and 
the O — C measurements were fixed to an appropriate location on the O — C diagram. 
Due to the large observational gap between 1994 and 1997, the measurements from 1994 
remained unconstrained. The appropriate location in the O — C diagram of the 1994 O — C 
measurements is further addressed in section 13.51 No unambiguous trend was apparent in 
the O — C diagrams of pulsation frequencies C, B, J, and K and they were removed from 
analysis. The O — C diagrams of pulsation frequencies A, D, E, G, H, and I are shown in 
Figures [T]|6l 



All six O — C diagrams (Figures [T]|n]) show some sort of oscillatory behavior. A sinusoid 
plus an offset, linear, and parabolic term were chosen as a suitable model. The model is 
mathematically described as 



Instead of fitting this nonlinear model with some initial guesses to t^, 6 fobs-, the rate 



the sinusoidal frequency variation 11, and the absolute phase of the sinusoidal frequency 
variation 0, a grid of linear fits can be calculated for a range of 11. By removing any 
dependence on initial guesses for the parameters, it is guaranteed that the absolute 
minimum in over the chosen range of 11 will be found. A Monte Carlo simulation 
was performed to calculate the statistical likelihood of the possible locations of the 1994 



3.5. Modeling the O — C Variations 





of period change P, the amplitude of the sinusoidal frequency variation a, the period of 
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150 



100 




-150' ^ ^ ^ ^ ' 

1995 2000 2005 2010 2015 

Time (Years) 

Fig. 1. — O — C of pulsation frequency A. The dashed hues are located at 0, P, and —P. 
The blue line is the best model fit for 11 = 12.9yrs. The red lines indicate the boundaries 
of the one sigma likelihood prediction of the model not including the 1994 data. Note that 
a first order polynomial has been removed from the data. 
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250 r 

200 - 
150 - 
100 




-150 - 

-200 -. 

-250 ' ^ ^ ^ ^ 

1995 2000 2005 2010 2015 

Time (Years) 

Fig. 2. — O — C of pulsation frequency D. The dashed hues are located at 0, P, and —P. 
The blue line is the best model fit for 11 = 12.9yrs. The red lines indicate the boundaries 
of the one sigma likelihood prediction of the model not including the 1994 data. Note that 
a first order polynomial has been removed from the data. 
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300 r 
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100 - 




-100 - 
-200 - 



-300 ' ^ ^ ^ ' ' 

1995 2000 2005 2010 2015 

Time (Years) 

Fig. 3. — O — C of pulsation frequency E. The dashed hues are located at 0, P, and —P. 
The blue line is the best model fit for 11 = 12.9yrs. The red lines indicate the boundaries 
of the one sigma likelihood prediction of the model not including the 1994 data. Note that 
a first order polynomial has been removed from the data. 
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300 




-300 ' ^ ^ ^ ' ' 

1995 2000 2005 2010 2015 

Time (Years) 

Fig. 4. — O — C of pulsation frequency G. The dashed hues are located at 0, P, and —P. 
The blue line is the best model fit for 11 = 12.9yrs. The red lines indicate the boundaries 
of the one sigma likelihood prediction of the model not including the 1994 data. Note that 
a first order polynomial has been removed from the data. 
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-200 - 



-300 ^ ^ ^ ^ ' ^ 

1995 2000 2005 2010 2015 

Time (Years) 



Fig. 5. — O — C of pulsation frequency H. The dashed hues are located at 0, P, and —P. 
The blue line is the best model fit for 11 = 12.9yrs. The red lines indicate the boundaries 
of the one sigma likelihood prediction of the model not including the 1994 data. Note that 
a first order polynomial has been removed from the data. 
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-400 ' ^ ^ ^ ' ' 

1995 2000 2005 2010 2015 

Time (Years) 



Fig. 6. — O — C of pulsation frequency I. The dashed hues are located at 0, P, — P, and 2P. 
The blue line is the best model fit for 11 = 12.9yrs. The red lines indicate the boundaries 
of the one sigma likelihood prediction of the model not including the 1994 data. Note that 
a first order polynomial has been removed from the data. 
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measurements given the model and the 1997-2012 measurements. The one sigma hkehhood 
boundaries are indicated in Figures [T]ini It was concluded that for all frequencies there 
was only one reasonable choice for the possible values of the 1994 O — C measurements 
and that all ambiguity could be considered resolved for the given model. The 1994 O — C 
measurements were then fixed to these values in each O — C diagram. Fits to Equation [3] 
were then calculated over a range of 11 for each of the O — C diagrams. A plot of reduced 
for the fit of each of the six O — C diagrams as a function of 11 is shown in Figure [3 Four 
of the six O — C diagrams have a strong periodicity at around 13 years while the other two 
are not particularly sensitive to the choice of period, but are consistent with a periodicity 
of 13 years. The minima of the total reduced for all six O — C diagrams is at 11 = 12.9 
years. This value of 11 was used to fit each O — C diagram. The best fits are overlayed in 
Figures [T] - E] and Table H] summarizes the resultant parameters. The large value of reduced 
for some of the fits is disturbing, but the ability of the model to reproduce the overall 
trend in the data is excellent. To further investigate, 11 was then allowed to vary for each 
individual O — C . This lead to a worse overall reduced x^, an indicator that the oscillatory 
behavior in each O — C is represented appropriately by the same 11. Fitting the data with 
high order polynomials did not result in a dramatic improvement in reduced x^. An 8th 
order polynomial, for example, actually increased reduced for two of the O — C diagrams 
and only led to an overall improvement of 7% over Equation [3l Equation [3] is dramatically 
more constrained than a polynomial with a similar number of degrees of freedom so this 
strongly supports the ability of Equation [3] to represent the overall trends in the data. 
The high values of reduced mean that either the error bars in the data are somehow 
underestimated or there are high frequency variations in O — C that are undersampled. To 
investigate, the Fourier transforms of the averaged residuals and absolute averaged residuals 
(see Figure [HI) were calculated. There was no statistically significant or notable power at 
any frequency, although the residuals appear correlated. It is concluded that there are 
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likely unresolved variations in phase, but these variations do not dilute the effectiveness of 
Equation [3] in modeling the overall variations. 

4. Discussion 

4.1. Planetary Detection with the Pulsation Timing Method 

The sinusoidal component of the variations shown in Figures [T]l6] do not share the 
same amplitude and phase, so are clearly not due to the effects of a planetary companion. 
Had only one of the pulsation frequencies in this star been analyzed, this paper could 
very well be announcing the first detection of a planet around a white dwarf. The 
discovery that multiple O — C diagrams can show similar periodic behavior when not in the 
presence of a planetary companion is alarming. Pulsation timing based white dwarf planet 
detection should now require supplemental confirmation, even when similar pulsation timing 
variations are observed in multiple pulsation frequencies. Until the physical mechanisms 
behind these variations are identified as specific to white dwarfs, these observations cast 
at least some shadow of doubt on the effectiveness of pulsation timing to reliably detect 
planetary companions to subdwarf B pulsators, and raise at least some skepticism as to 
the existence the subdwarf B planet V391 Pegasi b. We note that this does not affect the 
ability of the pulsation timing method to exclude the existence of planets around white 
dwarfs and other variables. The lack of known white dwarf planets is gaining statistical 
significance and further timing observations of pulsating white dwarfs could place strong 
constraints on post main sequence stellar evolution. 
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Table 3: Precision of frequency measurements from bootstrapping of data from September 
2003 to August 2004. Standard errors are quoted in parenthesis. 



Label 


Period [s] 


Frequency [uhz] 




A 


134 


7452. 249(. 001) 


32 


B 


195 


5128.598(.001) 


32 


C 


204 


4902.176(.0009) 


35 


D 


205 


4887.8491(.0009) 


35 


E 


257 


3893.2482(.0002) 


160 


G 


281 


3559.0017(.0003) 


106 


H 


287 


3489.056(.001) 


29 


I 


333 


2998.7144(.0007) 


45 


J 


350 


2852.439(.002) 


19 


K 


525 


1903.464(.001) 


24 



Table 4: Resultant fit parameters from fitting Equation [3] to the six O — C diagrams. Note 
the high values of reduced fit means that some of the formal quoted errors are unreliable, 
even though the model reproduces the overall trends in the data very well. 



Label 


Period [s] 


P [10-14] 


a [10-9i] 


$ [Deg] 




A 


134 


-28(1) 


4.0(.l) 


160(2) 


3.3 


D 


205 


-1(2) 


1.33(.09) 


187(4) 


6.6 


E 


257 


13.0(.8) 


2.48(.02) 


0(1) 


14 


G 


281 


-131(1) 


6.44(.03) 


166.6(.2) 


26 


H 


287 


88(6) 


6.8(.l) 


1(1) 


2.8 


I 


333 


313(4) 


7.31(.07) 


-1(1) 


13 
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Fig. 8. — Residuals of the fit of Equation |3] to the O — C diagrams of D,E,G,H and I. The 
black points are the weighted average residuals for each data set while the red points are the 
weighted average of the absolute value of the residuals. 
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4.2. The Combination Frequency A 



Sullivan et al 



(|2008[ ) identified the frequency of A to be an additive combination of 



the frequencies of E and G, a well known occurrence in white dwarf pulsators (see e.g. 



Provencal et al. 



20091 ). Any perturbation to the frequency of E or G should result in an 
identical perturbation to the frequency of A. The frequency perturbations found by fitting 
Equation [3] to the O — C diagram of E, G, and A can be compared, see Table [51 The 
agreement is outstanding. This is strong evidence that Equation [3] is effectively modeling 
the the frequency perturbations and that the parabolic and sinusoidal components of 
the perturbation are linearly independent. It is also the most conclusive evidence that 
combination fre quencies a re indeed exactly equal to the additive frequency of the parents 



as predicted by 



Wu 



mm . 



4.3. An Asteroseismologic Interpretation of the Results 



The frequency of these pulsations are dependent on many physical quantities, including 
the star's rotational velocity, magneti c field strength, and magnetic field geometry (for an 



in depth review, see 



Unno et al. 



19791 ). Any redistribution of angular momentum or change 



in magnetic field strength (or geometry) would change the observed frequencies. The sign of 



Table 5: Comparison of the modeled variations of the combination frequency A to its parents 
E and G. 



Parameter 


A 


E + G 


/ = -^[^] 
«[^] 
(p [deg] 


14.6(.2) 
4.0(.l) 
160(2) 


15.3(.5) 
4.07(.03) 
158.7(.2) 
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the first order frequency perturbation caused by changes in rotation is well known to depend 
solely on whether the pulsations are moving with or against rotation. In other words, if 
one pulsation is moving with rotation and another against it, a redistribution of angular 
momentum would perturb the frequency of the two pulsations in the opposite direction. 
The variations in pulsation frequencies E, H, and I have positive P and degrees while 
pulsation frequencies D and G have negative P and close to 180 degrees. We speculate 
that one of these groups of pulsation frequencies could be moving with rotation and the 
other against rotation and that the observed variations in O — C may be caused by changes 
in the star's rotation profile. Additionally, the long period g-mode pulsation frequencies in 
white dwarfs are well known to be more affected by physical properties at the surface of the 
star than short period g-mode pulsation frequencies. The observed frequency perturbations 
appear to affect the longer period pulsations more, a hint that the corresponding physical 
processes responsible for these frequency variations may be near the surface of the star. 
We again emphasize the speculative nature of these interpretations. The possible effects of 
magnetic field variations and quantitative modeling of these frequency variations will be 
investigated in the primary author's doctoral dissertation. 
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